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1. Introduction 

The newly created ultracold trapped Fermi gases with tunable atomic scattering length 
[1-20] in the vicinity of a Feshbach resonance offer the possibility to study highly corre- 
lated many-body systems including the cross-over from the Bardeen-Cooper-Schrieffer (BCS) 
phase to the Bose-Einstein condensate (BEC) of molecules. Various investigations based on 
the hydrodynamic theory have appeared recently [21-31]. 

The purpose of this letter is to construct a simple, semi-analytical and nearly exact 
formulae for hydrodynamic frequencies. Since the collective frequencies can be measured 
with high precision, these formulae will provide a simple quantitative tool for the analysis 
of experimental data in the hydrodynamic regime. 

2. Hydrodynamic theory 

Our starting point is the quantum hydrodynamic theory [21-23] for a dilute two-component 
Fermi gas in a trap potential V ext (r) = (m/2)(iu'j_(x 2 + y 2 ) + uj 2 z 2 ) 

ih— = -— V 2 ^ + V ext V + V XC V, (1) 
ot 2m 

where 

V (r A - [ d( ne ( n )) l 

VxcV i^) — [ 7^ \n=n{r,t)i 

e is the ground state energy per particle of the homogeneous system and n is the density, 
n(r, t) =| \l/(r, t) | 2 , normalized to the total number of atoms, / n(f,t)d 3 r = N. It is useful 
to rewrite Eq.(l) in a form 

% + V(nv) = 0, (2) 
dv d(ne(n)) 1 , h 2 1 „, ^ N 

where v is the velocity field, which for ^ = e l ^ r ' t ^n 1 ^ 2 (f,t) can be written as v = (h/m)V(f). 

It can be proved [21] that every solution of the Eqs.(2,3) is a stationary point correspond- 
ing to the Lagrangian density 

h 2 - h 2 

C = Hn + ^(Vv 7 ^) 2 + ^ n W) 2 + e ( n ) n + V ^ n - ( 4 ) 

It was shown in Refs. [21-23] that for experimental conditions of Refs. [14,16,17] the quantum 
pressure term in Eqs.(3,4) can be neglected. For the reminder of this Letter we will use this 
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hydrodynamic approximation. For the harmonic trap a trial function in the scaling ansatz 
is taken as [21-23,26,28] 

(f>(f,t) = Mt) + (V(2fi)EA(fX ? > 

i=i 

n(f,t) = no(xi/bi(t))/K(t), 

where n(t) = Ylj bj and the Hamilton principle, 5 J dt J Cod 3 r = 0, gives the following equa- 
tions for the scaling parameters b{ [21-22] 

bi J[n 2 rfe(n)/dn]„ = „ o(TO dV 

Expanding Eqs.(5) around equilibrium (pi = 1) leads to the following result for the M = 
modes frequencies, to^ in the scaling approximation 

"S? = ^[Vs ± v/^-SA^SC + S)] 1 / 2 , (6) 

where r) s = 4 + 2( s + 3A 2 + ( S A 2 , ( s — J nf ) d 2 e/(dnl)d 3 r/ J n\dej (dn )d 3 r, A = cu z /cu±, and ± 
refer to the transverse and axial mode, respectively. 

The hydrodynamic equations after linearization take the form 



where 5n(r, t) is the change in the density profile with respect to the equilibrium configura- 
tion. If we consider oscillations with time dependence 5n oc exp(icut), Eq.(7) can be reduced 
to a Hermitian equation [27] 



where [/ >= d ( re °^"°^ \Sn >, L = — — Vn V and the equilibrium density, n , is given by 

ClTIq 771 

equation 

dn 

where \i is the chemical potential, in the region where no(r) is positive and no(r) = outside 
this region. 

We note here that d 2 (n e(n ))/dn 2 ) is positive, since the sound velocity for the homoge- 
neous case is given by c 2 = (no/m)d 2 (noe(no)) / dn 2 ^. 
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3. Equation of state 

For the negative S-wave scattering length between the two fermionic species, a < 0, in 
the low-density regime, fcf | a |< 1, the ground state energy per particle , e(n), is well 
represented by an expansion in power of k F \ a | [32] 

e(n) = 2E F [^- - ^-k F | a | +0.055661(A; F | a |) 2 - 0.00914(A; F | a |) 3 + ...], (10) 

10 07T 

where E F = h 2 k F /(2m) and k F = (37r 2 n) 1//3 . In the opposite regime, a — > — oo (the Bertsch 
many-body problem, quoted in Refs.[33]), e(n) is proportional to that of the non- interacting 
Fermi gas 

e(n) = (l+(3)^^, (11) 
10 m 

where a universal parameter f3 [10] is estimated to be (3 = —0.56 [34]. The universal limit 
[10,34-37] is valid at least in the case where the width of the Feshbach resonance is large 
compared to the Fermi energy as in the cases of 6 Li and 40 K. 

In the a — > +0 limit the system reduces to the dilute Bose gas of dimers 

e(n) = E F {-l/{k F af + a m k F /(Qn) + ...,) (12) 

where a m is the boson-boson scattering length, a m ~ 0.6a [38]. 

A simple interpolation of the form e(n) ~ E F P(k F a) with a smooth function P(x) was 
considered in several papers. In Ref. [21] a [2/2] Pade approximant has been proposed for 
the function P(x) for the the case of negative a 

P( X) = 3 -- 2 , (13) 

where ^ = 0.106103, 5 2 = 0.187515, 5 3 = 2.29188, 5 A = 1.11616. Eq.(13) is constructed 
to reproduce the first four terms of the expansion (10) in the low-density regime and also 
to exactly reproduce results of the recent Monte Carlo calculations [34], (3 = —0.56, in the 
unitary limit, k F a — > — oo. 

For the positive a case ( the interaction is strong enough to form bound molecules with 
energy E mo {) we have considered in Ref. [22] a [2/2] Pade approximant 

P( T \ - E ™i | aiX + ^ n A \ 
1 >~ 2E F + l + a 3 x + a 4 x^ [ } 

where parameters ctj are fixed by two continuity conditions at large x, 1/x — > 0, and by 

two continuity conditions at small x, oti = 0.0316621, a 2 = 0.0111816, a 3 = 0.200149, and 

a A = 0.0423545. 
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In Ref. [39] a Pade approximation has been considered for the chemical potential. Authors 
of Ref. [28] have used a model for P(x) , interpolating the Monte Carlo results of Ref. [36] across 
the unitary limit and limiting behaviors for small \x\. We note here also the BCS mean-field 
calculations of Ref. [28]. 

4. Schwinger variational principle 

The Schwinger variational principle (SVP) [40-43] can be generalized to the case of Eq.(8). 

Since the entire treatment is based on the equivalence of the SVP and the method of 
separable representation, we briefly describe this method. Let us consider the symbolic 
identity 

L = LL- l L = Y, L \i >< il^lj >< ( 15 ) 

where \i > is a complete set. Truncating the summation over the complete set we obtain a 
separable approximation 

L<«> = Y.L\ Xl >d- 1 <Xj \L, (16) 

where d^ =< Xi\L\Xj >■ 

We note that Eq.(16) represents an interpolation process, since L^\xi >= L\Xi > an d 
< Xj \L^ =<Xj\L. 

Substituting from Eq.(16) into Eq.(8), we obtain 

u2 fMn )) nf >= L(9)|/ >= f L{Xt > rf ! K > 

We seek a solution of Eq.(17) in the form 

* rf 2 (n e(n )) 

\f>=2^ °i 3T2 >» ( 18 ) 



then q's are defined from equations 



J2B ik (u 2 )c k = 0, (19) 
k=i 



B lk {u 2 ) =< Xi\(" 2 L - L K Z2 , L )\Xk > (20) 



where 

d 2 ( n e(n )) 
dn\ 

and frequencies cu are determined from the condition of vanishing of the determinant of the 
matrix Bik{oJ 2 ). 

detB ik {u 2 ) = (21) 



Defining the Schwinger functional Isvp[x] by 

| L d2( (no)) 

= ' < ^ x> ' ■ («) 

we get isvp[/] = w 2 , where / is the solution of Eq.(8). Introducing the function X — 
J2r=i c rXr, where ci,C2,...c g are g variable parameters, we see that the functional Isvp is 
stationary if dlsvp/dc r = 0. The later equations coincide with Eqs.(19), and therefore the 
approximate solution of the Eq.(8), which is based on the SVP, is equivalent to the exact 
solution of the Eq.(8) with separable lM\ 

Operator L is clearly positive, which means that < u\L\u >> for all u, but is not 
positive definite, since < u\L\u >= for some u ^ 0. It can easily be seen that 

< u\(L-L {q) )\u >> 0, (23) 

for all u. Indeed, 

g q 

J =< (u + J2c iX i)\L\(u + J2c jXj ) >> 
i=i j=i 

for all Cj. We choose the q's from the conditions dJ/dci = 0, and then obtain 

J=< u\(L-L {q) )\u >> 

for all q. 

Since the problem is solved by replacing L by and since the operator (L — L^) is 
positive, the SVP leads to the approximate lower bounds for u up to the second order of 
\L-Lti>\. 

For the most interested case of M = modes, we can put in Eq.(16) q — 2, Xi = ( a;2 +l/ 2 ) 
and X2 = z 2 , which give 

"1 VP = ^[vsvp ± ^ 2 svp-m%svP-l)] l/ \ (24) 

where r) SV p = 6( S vp+\ 2 (3( svp +1), ( SV p = - J nffi^^dx/ f n x 4 dx, x = y/x' 2 + y' 2 + X 2 z 2 
and ± signs refer to the transverse and axial mode, respectively. 

It is easy to show that Eq. (24) gives exact solutions for frequencies of the breathing 
modes for the polytropic equation of state, e(n) ~ ri 1 . In Ref. [24], on the basis of a gener- 
alization of the Hylleraas-Undheim method, we have constructed rigorous upper bounds to 
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the collective frequencies for the radial and the axial breathing mode of the Fermi gas under 
harmonic confinement in the framework of the hydrodynamic theory 

^VV _ Wjupper i \JVupper ~ ^^Cupper {^Cupper ~ 9)] ^ , (25) 

Y ^Supper " 

where T] upper = (3 + 4:X 2 )( upper - (3 + 6A 2 ), C« PP er = hh/Ih and ^ = / x l n (x)dx. 

We expect that the difference between u;| yp and u ^ pper characterizes the error and is not 
very sensitive to the functional form of e(n). From Table 1 one can see that this difference is 
order of 10~ 5 . We note that the comparison with the scaling approximation, Eq.(6), shows 
that the absolute precision of the scaling approximation is about 1CT 3 that agrees with 
Ref.[31]. 

In Fig.l, we have compared the hydrodynamic predictions for oj + with experimental data 
[17]. There is a very good agreement with experimental data [17] near the unitary limit. We 
note here that two experimental results [17] and [16] (not shown in Fig.l) for oo + are still 
about 10% in disagreement with each other, which is not fully understood yet. 

To calculate £'s, we have used the very fast converged expansion of Ref.[24] 

n (f) « (1 - (3V ext (^)) 1/{2 - p) EQ^UOr, 

i=0 

where parameters ft, p and q are fixed by requiring that n (r) must satisfy a variational 
principle 5 J n (V ext + e(n ))d 3 r = with a subsidiary condition / n d 3 r = N. 

4. Summary 

We have generalized the Schwinger variational method for the trapped strongly interact- 
ing atoms in hydrodynamic regime and we have constructed semi-analytical and extremely 
accurate formulae for hydrodynamic collective frequencies. These formulae are very useful 
since they provide an easy and simple quantitive tool for the analysis of experimental data 
for trapped condensed gases without relying on complex and extensive computations. 
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Table 1. The transverse and axial frequencies in units of uo± and u z , respectively, in 
the BCS region as a function of the dimensional parameter X = (N l l & a/ a ho )~ l . The trap 
parameter A is assumed to be 0.045613. The [2/2] Pade approximation of Refs. [21,22] is used 
for the energy per particle e(n). 



X 


upper 


. ,SVP 


^upper 


co 3 ^ 




UJ S _ 


-0.1 


1.8160 


1.8160 


1.5470 


1.5470 


1.8193 


1.5477 


-0.3 


1.8015 


1.8015 


1.5438 


1.5438 


1.8082 


1.5453 


-0.5 


1.7931 


1.7931 


1.5419 


1.5419 


1.8002 


1.5435 


-0.7 


1.7886 


1.7886 


1.5409 


1.5409 


1.7947 


1.5423 


-0.9 


1.7867 


1.7867 


1.5405 


1.5405 


1.7910 


1.5414 


-1.1 


1.7861 


1.7861 


1.5403 


1.5403 


1.7886 


1.5409 


-1.3 


1.7865 


1.7865 


1.5404 


1.5404 


1.7871 


1.5406 


-1.5 


1.7873 


1.7873 


1.5406 


1.5406 


1.7863 


1.5404 


-1.7 


1.7884 


1.7884 


1.5409 


1.5409 


1.7860 


1.5403 


-2.0 


1.7902 


1.7902 


1.5413 


1.5413 


1.7861 


1.5403 
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Fig. 1. Radial breathing mode frequency uj + in the BCS region as a function of the 
dimensional parameter X = (iV^a/a^o) -1 (the solid line). The dashed line represent the 
scaling approximation, Eq.(6). The circular dots with error bars are the experimental results 
given by the Duke University group [17]. Everything is measured in units of u±. 
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